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ABSTRACT 

We have studied the relationships between the electronic structure and the antiviral activity againat H1N1, 
enterovirus 71 and Coxsackie B3 pathogens. The local atomic reactivity indices were obtained at the B3LYP/6-31G (d,p) 
level. For two of them (coxsackievirus B3 and H1N1) we obtained statistically significant equations that led to the 
correponding pharmacophores displaying atomic sites that can be sustituted to obtain molecules endowed with higher 
inhibitory activity. The results for enterovirus 71 strongly suggest that the molecules seem to have more than one 
mechanism of action. 
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INTRODUCTION 

Influenza A (H1N1) virus is the subtype of influenza A virus that was the most common cause of human influenza 
in 2009 and is associated with the 1918 outbreak that killed 50 to 100 million people. Enterovirus 71 (EV71) is a virus 
notable for its etiological role in epidemics of severe neurological maladies in children. EV71 occasionally causes 
polio-like syndrome permanent paralysis. Coxsackie B3 has been found to be one of the chief causes of some weakening or 
life-threatening maladies, such as viral myocarditis. Several groups of molecules acting again these pathogens have been 
synthezised and tested (see for example [1-13]). 

One of the dreams of a quantum pharmacologist is the opportunity to study a set of molecules having more than 
one biological action or having a biological effect on different systems. Recently, a group of of pyrazine-l,3-thiazine 
hybrid analogues was synthezised and tested against influenza A (H1N1) virus, enterovirus 71 (EV71), and coxsackievirus 
B3 (CVB3) [14]. Given the importance of finding more molecular agents able to combat these pathogens, we presente here 
the results of the application of a formal method relating the electronic structure with the biological activity for a group of 
molecules displaying activity against the three abovementioned viruses, of the paper should explain the nature of the 
problem, previous work, purpose, and the contribution of the paper. The contents of each section may be provided to 
understand easily about the paper. 
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METHODS, MODELS AND CALCULATIONS 

As the method employed here has been widely discussed in many publications, we refer the reader to the literature 
[15-21]. Its application to molecule-site binding constants and to other biological activities has produced excellent results 
[22-49] (and references therein). Therefore, we shall discuss here only the results obtaines from its application. 
The selected molecules are shown in Figure 1 and Table 1. The biological activities analyzed here (Table 1) are the 
inhibition of the virus-induced cytopathic effect (CPE) in Madin-Darby canine kidney cells infected with the H1N1 virus 
(expressed as IC 50 ), and the inhibitory activity against Enterovirus 71 (EV71), expressed as IC 50 ) and coxsackievirus B3 
(CVB3, expressed as IC 50 ). 



Figure 1: Pyrazine-1, 3-Thiazine Hybrid Analogues 


Table 1: Pyrazine-1, 3-Thiazine Hybrid Analogues and Biological Activities 


Mol. 

Mol. 


r 3 

r 4 

log(IC 5 ») 

H1N1 

log(IC 50 ) 

EV71 

log(ICso) 

CVB3 

1 

3a 

H 

H 

H 

1.52 

1.19 

1.29 

2 

3b 

H 

H 

OMe 

1.05 

0.73 

1.15 

3 

3c 

H 

OMe 

H 

0.96 

0.96 

1.05 

4 

3d 

OMe 

H 

H 

0.73 

1.25 

1.48 

5 

3e 

H 

H 

Me 

1.31 

1.05 

1.28 

6 

3f 

H 

Me 

H 

1.23 

0.95 

0.92 

7 

3g 

Me 

H 

H 

1.27 

0.35 

0.80 

8 

3h 

H 

H 

OH 

1.38 

1.26 

1.46 

9 

3i 

H 

OH 

H 

1.47 

1.21 

1.44 

10 

3j 

OH 

H 

H 

1.50 

0.39 

1.34 

11 

3k 

H 

H 

no 2 

1.60 

1.55 

1.26 

12 

31 

H 

no 2 

H 

1.65 

1.59 

1.28 

13 

3m 

no 2 

H 

H 

1.62 

0.82 

0.95 

14 

3n 

H 

H 

Cl 

1.73 

1.55 

1.48 

15 

3o 

H 

Cl 

H 

1.76 

1.69 

1.59 

16 

3p 

Cl 

H 

H 

1.73 

1.74 

1.51 

17 

3q 

H 

H 

F 

1.94 

1.29 

1.53 

18 

3r 

H 

F 

H 

1.97 

1.20 

1.37 

19 

3s 

F 

H 

H 


1.38 

1.62 


The electronic structure of all molecules was calculated within the Density Functional Theory (DFT) at the 
B3LYP/6-31g(d,p) level with full geometry optimization [50]. The Gaussian suite of programs was used [51]. All the 
information needed to calculate numerical values for the local atomic reactivity indices was obtained from the Gaussian 
results with the D-Cent-QSAR software [52]. All the electron populations smaller than or equal to 0.01 e were considered 
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as zero [19]. Negative electron populations coming from Mulliken Population Analysis were corrected as usual [53]. 
Orientational parameters were taken from the literature [54, 55]. Since the resolution of the system of linear equations is 
not possible because we have not enough molecules, we made use of Linear Multiple Regression Analysis (LMRA) 
techniques to find the best solution. For each case, a matrix containing the dependent variable (the log(IC 50 ) of each case) 
and the local atomic reactivity indices of all atoms of the common skeleton (see below) as independent variables was built. 
The Statistica software was used for LMRA [56]. We worked with the common skeleton hypothesis stating that there is a 
definite collection of atoms, common to all molecules analyzed, that accounts for nearly all the biological activity [18]. 
The action of the substituents consists in modifying the electronic structure of the common skeleton and influencing the 
right alignment of the drug throughout the orientational parameters. It is hypothesized that different parts or this common 
skeleton accounts for almost all the interactions leading to the expression of a given biological activity. The common 
skeleton for this case is shown in Figure 2. 



Figure 2: Common Skeleton Numbering 


RESULTS 

Results for the Inhibition of the Virus-Induced Cytopathic Effect in Madin-Darby Canine Kidney Cells Infected 
with the H1N1 Virus 

The best equation obtained is: 

log(IC 50 )=2.07+0.09S* (LUMO+2)*-l ,39F 3 (HOMO-2)*+l . 1 3F 16 (LUMO+2)*- 
-1.07F 18 (HOMO-2)*+0.81F 10 (LUMO+2)* 

With n=16, R=0.98, R 2 =0.97, adj-R 2 =0.95, F (5, 10) =58.53 (p<0.000001) and a Std. error of estimate of 0.07. 
No outliers were detected and no residuals fall outside the +2cr limits. Here, Sn N (LUMO+2)* is the nucleophilic 
superdelocalizability of the third lowest empty MO localized on atom 11, F 3 (HOMO-2)* is the Fukui index of the third 
highest occupied MO localized on atom 11, Fi 6 (LUMO+2)* is the Fukui index of the third lowest empty MO localized on 
atom 16, Fi 8 (HOMO-2)* is the Fukui index of the third highest occupied MO localized on atom 18 and Fi 0 (LUMO+2)* is 
the Fukui index of the third lowest empty MO localized on atom 10. Tables 2 and 3 show the beta coefficients, the results 
of the t-test for significance of coefficients and the matrix of squared correlation coefficients for the variables of Eq. 1 . 
There are no significant internal correlations between independent variables (Table 3). Figure 3 displays the plot of 
observed vs. calculated log (IC50). 
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Table 2: Beta Coefficients and t-Test for Significance of Coefficients in Eq. 1 


Variable 

Beta 

t(10) 

p-level 

S,i (LUMO+2)* 

0.88 

10.95 

<0.000001 

F 3 (HOMO-2)* 

-0.37 

-5.38 

<0.0003 

F 16 (LUMO+2)* 

0.56 

7.86 

<0.00001 

F lg (HOMO-2)* 

-0.27 

-3.94 

<0.003 

F 10 (LUMO+2)* 

0.16 

2.74 

<0.02 


Table 3: Matrix of Squared Correlation Coefficients for the Variables in Eq. 1 



SxALUMO+2)* 

F 3 (HOMO-2)* 

F 16 (LUMO+2)* 

F 18 (HOMO-2)* 

F 3 (HOMO-2)* 

0.20 

1 



F 16 (LUMO+2)* 

0.22 

0.05 

1 


F 18 (HOMO-2)* 

0.03 

0.05 

0.06 

1 

F 10 (LUMO+2)* 

0.01 

0.02 

0.02 

0.01 



Figure 3: Plot of Predicted Vs. Observed Log (IC 50 ) Values (Eq. 1). 

Dashed Lines Denote the 95% Confidence Interval 

The associated statistical parameters of Eq. 1 indicate that this equation is statistically significant and that the 
variation of the numerical values of a group of five local atomic reactivity indices of atoms of the common skeleton 
explains about 95% of the variation of log (IC 50 ). Figure 3, spanning about 1.4 orders of magnitude, shows that there is a 
good correlation of observed versus calculated values. 

Results for the Inhibition of Enterovirus 71 (EV71) 

The best equation obtained is: 

log(IC 50 )=-49.83+0.23S^ 4 (LUMO+2)*-0.006S, N ,(LUMO+2)*-84.71Q 19 ... 

( 2 ) 

With n=19, R=0.96, R 2 =0.92, Adj-R 2 =0.91, F (3, 15) =59.932 (p<0.000001) and a Std. error of estimate of 0.12. 
No outliers were detected and no residuals fall outside the +2a limits. Here, S 2 4 N (LUMO+2)* is the nucleophilic 
superdelocalizability of the third lowest empty MO localized on atom 24, Si 3 N (LUMO+2)* is the nucleophilic 
superdelocalizability of the third lowest empty MO localized on atom 13 and Q 19 is the net charge of atom 19. Tables 4 and 
5 show the beta coefficients, the results of the t-test for significance of coefficients and the matrix of squared correlation 
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coefficients for the variables of Eq. 2. There are no significant internal correlations between independent variables 
(Table 7). Figure 4 displays the plot of observed vs. calculated log (IC 50 ). 

Table 4: Beta Coefficients and t-Test for Significance of Coefficients in Eq. 2 


Variable 

Beta 

t(15) 

p-level 

S 24 n (LUMO+2) * 

0.73 

9.78 

< 0.000000 

S 13 N (LUMO+2)* 

-0.56 

-7.30 

<0.000003 

Ql9 

-0.25 

-3.22 

<0.006 


Table 5: Matrix of Squared Correlation Coefficients for the Variables in Eq. 2 



S 24 N (LUMO+2)* 

S 13 N (LUMO+2)* 

S 13 (LUMO+2)* 

0.05 

1 

Q 19 

0.05 

0.11 



Figure 4: Plot of Predicted Vs. Observed Log (IC 50 ) Values (Eq. 2). 

Dashed Lines Denote the 95% Confidence Interval 

The associated statistical parameters of Eq. 2 indicate that this equation is statistically significant and that the 
variation of the numerical values of a group of three local atomic reactivity indices of atoms of the common skeleton 
explains about 91% of the variation of log (IC 50 ). Figure 4, spanning about 1.3 orders of magnitude, shows that there is a 
good correlation of observed versus calculated values. 

Results for the Inhibition of Coxsackievirus B3 (CVB3) 

The best equation obtained is: 

log(IC 50 )=57.98+0.40Sf 8 (HOMO-2)*-1.10F 17 (LUMO+2)*-0.10S2 4 (LUMO+l)*- 

-4-84Q 16 +125.9 1Q 2 +0.28Q 14 (3) 

With n=19, R=0.97 R 2 =0.94, adj-R 2 =0.91, F(6,12)=29.649 (p<0.00001) and a Std. error of estimate of 0.07. No 
outliers were detected and no residuals fall outside the ±2o limits. Here, Sig E (HOMO-2)* is the electrophilic 
superdelocalizability of the third highest occupied MO localized on atom 18, F| 7 (LUMO+2)* is the Fukui index of the 
third lowest empty MO localized on atom 17, S 24 N (LUMO+1)* is the nucleophilic superdelocalizability of the second 
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lowest MO localized on atom 24, Q 16 is the net charge of atom 16, Q 2 is the net charge of atom 2 and Q 14 is the net charge 
of atom 14. Tables 6 and 7 show the beta coefficients, the results of the t-test for significance of coefficients and the matrix 
of squared correlation coefficients for the variables of Eq. 3. There are no significant internal correlations between 
independent variables (Table 7). Figure 5 displays the plot of observed vs. calculated log (IC 50 ). 

Table 6: Beta Coefficients and t-Test for Significance of Coefficients in Eq. 3 


Variable 

Beta 

t(12) 

p-level 

S 18 k (HOMO-2)* 

0.55 

6.16 

<0.00004 

F 17 (LUMO+2)* 

-0.63 

-6.85 

<0.00002 

S 24 n (LUMO+1)* 

-0.46 

-5.46 

<0.0001 

Q 16 

-0.36 

-3.40 

<0.005 

Q 2 

0.40 

3.66 

<0.003 

Q 14 

0.21 

2.23 

<0.05 


Table 7: Matrix of Squared Correlation Coefficients for the Variables in Eq. 3 



Si 8 e (HOMO-2)* 

F 17 (LUMO+2)* 

S 24 n (LUMO+1)* 

Qig 

O 2 

F 17 (LUMO+2)* 

0.04 

1 




S 24 n (LUMO+1)* 

0.04 

0.11 

1 



Q 16 

0.14 

0.31 

0.23 

1 


q 2 

0.20 

0.12 

0.13 

0.31 

1 

Q 14 

0.00 

0.00 

0.08 

0.10 

0.25 



Figure 5: Plot of Predicted Vs. Observed Log (IC 50 ) Values (Eq. 3). 
Dashed Lines Denote the 95% Confidence Interval 


The associated statistical parameters of Eq. 3 indicate that this equation is statistically significant and that the 
variation of the numerical values of a group of six local atomic reactivity indices of atoms of the common skeleton 
explains about 91% of the variation of log (IC 50 ). Figure 5, spanning about 0.8 orders of magnitude, shows that there is a 
good correlation of observed versus calculated values. 
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LOCAL MOLECULAR ORBITALS 

Tables 8 to 10 show the local MO structure of atoms 2-4, 10-13, 16-18, 20, 23 and 24 (see Figure 2). 
Nomenclature: Molecule (HOMO) / (HOMO-2)* (HOMO-1)* (HOMO)* - (LUMO)* (LUMO+1)* (LUMO+2)*. 

Table 8: Local Molecular Orbitals of Atoms 2, 3, 4 and 10 


Mol. 

Atom 2 

Atom 3 

Atom 4 

Atom 10 

1(82) 

80o81o82o- 

83ji847i857i 

79c80o81o- 

837i85jt877i 

80o81g82g- 

837i847i877t 

807i8l7i827t- 

83c847i867i 

2(90) 

85ji867i88o- 

91ji927i937i 

857i867i88o- 

9l7i927i937i 

88o89g90g- 

9l7i927i957i 

887i897i907i- 

9l7i927i947t 

3(90) 

867i87o88o- 

91ji927i937i 

867i87o88o- 

9l7i927i937i 

88o89g90g- 

9l7i927i957i 

88o897t907i- 

9l7i927i947i 

4(90) 

867i87o88o- 

91ji927i937i 

867i87c88a- 

9l7i927i937i 

87o88g90g- 

9l7i927i937t 

88a897t907i- 

9l7i927i937t 

5(86) 

84a85o86o- 

877i887i897i 

83o84c85o- 

877i88ji897i 

84g85g86g- 

877i887i9l7t 

84g85g86ti- 

877i887i907i 

6(86) 

83a84o86o- 

87jt887i897i 

827i83g84c- 

877i887i897i 

84g85g86g- 

877i887i9l7t 

847i857i867i- 

877i887i907i 

7(86) 

83a84a86o- 

877i887i897i 

827i83c84c- 
877i887i 8971 

84g85g86g- 
877i887i 8971 

847i857i867i- 
877188?! 8971 

8(86) 

81tc827i84o- 

877i887i897i 

817i82tt84o- 

877i887i897i 

84g85g86g- 

877i887i9l7t 

847i857i867i- 

877i887i9l7t 

9(86) 

827i83a84a- 

877i887i897i 

827i83o84o- 

877i887i897i 

84g85g86g- 

877i887i9l7t 

847i857i867i- 

877i887i907i 

10(86) 

83a84a86c- 

87jt887i897i 

827i83o84o- 

877i887i897i 

83g84g86g- 

877i887i897t 

847i85 718671 - 
877i88ti89o 

11(93) 

91o92o93o- 

95ji967i977i 

907i91o92o- 

957i967i977i 

91g92g93g- 

957i967il007i 

907i9l7i937t- 

957i967i977t 

12(93) 

91o92a93a- 

957i967i977i 

907i91g92o- 

957i967i977t 

91g92g93g- 

957i967i977t 

907i9l7i937t- 

957i967i977t 

13(93) 

91ji92o93o- 

95jr967i977i 

907i91ti92o- 

957i967i977i 

91ti92g93g- 

957i967i977i 

9l7i927i93o- 

957i967i977t 

14(90) 

88a89a90a- 

9l7i927i937i 

877i88g89g- 

9l7i937i957t 

88g89g90g- 

9l7t927i937c 

887i897i907i- 

9l7i927i947t 

15(90) 

88a89a90a- 

9l7r927i937i 

87g88g89g- 

9l7i937i957i 

88g89g90g- 

9l7i927i937i 

887i897i907i- 

9l7i927i947i 

16(90) 

88a89a90a- 

9l7i927i937i 

877i88g89g- 

9l7i927i937t 

88g89g90g- 

9l7i927i967i 

877i897i907i- 

9l7i927i947t 

17(86) 

84a85a86a- 

87ji887i897i 

837i84o85c- 

877i897i9l7t 

84g85g86g- 

877i887i9l7i 

84g857i86ti- 

877i887t9l7t 

18(86) 

84a85a86a- 

87ji887i897i 

83o84g85g- 

877i897i907i 

84g85g86g- 
877r887r91 7i 

84o85ct867i- 

877r887i907i 

19(86) 

84o85g86g- 

87ji887i897i 

837i84o85g- 

877i887i897i 

84o85g86g- 

877i887i927i 

847i857i867i- 

877i887i907i 
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Table 9: Local Molecular Orbitals of Atoms 11, 12, 13 and 16 


Mol. 

Atom 11 

Atom 12 

Atom 13 

Atom 16 

1(82) 

76a77o82Ti- 
8371 84 ti857i 

77a79a827t- 
8371 847i857i 

8 O 718 I 718271 - 
8371 84ti857i 

807i8l7t827i- 

86ti907i917i 

2(90) 

887i897i907i- 

9l7i927i937i 

887i897i907i- 

9l7i927i937i 

877i897i907i- 

927i937i947i 

887i897i907i- 

947i987i997i 

3(90) 

867i897i907i- 

9l7i927i937i 

877i897i907i- 

9l7i927i937i 

877i88o907i- 

917i92ti937i 

877i88a897i- 

937i947i957t 

4(90) 

887i897i907i- 

9l7i927i937i 

887i897i907i- 

91tt927i937i 

887i897i907i- 

927i937i947i 

887i897i907i- 

937i957i967i 

5(86) 

827i837i86ti- 

877i887i897i 

837i85o867i- 

877i887i897i 

837i857i86ti- 

877i887i897i 

847i857i86ti- 

907i947i957i 

6(86) 

81a827i867i- 

877i887i897i 

837i857i867i- 

877i887i897i 

847i857i86tt- 

877i887i897i 

837i84ti857i- 

897i907i947i 

7(86) 

827i837i86ti- 

877i887i897i 

837i85o867i- 

877i887i897i 

847i857i86ti- 

887i897i907i 

847i857i86ti- 

897i907i91tt 

8(86) 

847i857i86tt- 

877i887i897i 

847i857i86ti- 

877i887i897i 

847i857i86ti- 

887i897i907i 

847i857i86ti- 

907i91ti957i 

9(86) 

827i857i86ti- 

877i887i897i 

837i857i86ti- 

877i887i897i 

837i847i86ti- 

877i887t897i 

837i847i857i- 

897i907i91ti 

10(86) 

847i857i867i- 

877i887i897i 

847i857i867i- 

877i887i897i 

847t857i867i- 

887i897i907i 

847i857i867i- 

897i907i91ti 

11(93) 

917i92ti937i- 

947i957i967i 

90o927i937i- 

947i957i967i 

9l7i927i937i- 

947i977i987i 

897i907i91tc- 

947i957i967i 

12(93) 

917i92ti937i- 

947i957i967i 

90o927i937i- 

947i957i967i 

9l7i927i937i- 

947i957i967i 

897i907i91tc- 

947i957i967i 

13(93) 

917i92ti937i- 

947i957i967i 

90o927i937i- 

947i957i967i 

897i907i9l7i- 

947i957i967i 

887i897i907t- 

947i987i997i 

14(90) 

85o877t907i- 

91tt927i937i 

86o887i907i- 

91ti927i937i 

887i897i907i- 

917i92ti937i 

887i897i907i- 

937i947i957i 

15(90) 

85o887i907i- 

91ti927i937i 

85o887i907i- 

91ti927i937i 

887i897i907i- 

917i92ti937i 

877i887t897i- 

937i947i97o 

16(90) 

86o877i907i- 

91ti927i937i 

86o887i907i- 

91tt927i937i 

877t897i907i- 

917i92ti937i 

887i897i907i- 

92ti937i947i 

17(86) 

81a837i867i- 

877i887i897i 

82o847i867i- 

877i887i897i 

847i857i867i- 

877i887i897i 

847i857i867i- 

907i9l7i947i 

18(86) 

81o847i867i- 

877i887i897i 

837i847i867i- 

877i887i897i 

847i857i867i- 

877i887i897i 

837i847i857i- 

897i907i91ti 

19(86) 

82o837i867i- 

877i887i897i 

82o847i867i- 

877i887i897i 

837i857i867i- 

877i887i897i 

847i857i867i- 

897i907i91ti 


Table 10: Local Molecular Orbitals of Atoms 17, 18, 20 and 24 



Atom 17 

Atom 18 

Atom 20 

Atom 23 

Atom 24 

1(82) 

797i807i827i- 

847i857i867i 

807i8l7i827i- 

83tt847i857i 

78o79o80o- 

83o95o96o 

50o52o53o- 

87o88o89a 

63o64o77o- 

89c90o91o 

2(90) 

877i897i907i- 

927i937i947t 

877i897i907i- 

91tt927i937i 

84o86o88o- 

91ol04ol05a 

55a56o59o- 

95a96a97a 

61o70o85o- 

97o98o99o 

3(90) 

83897i907i- 

927i937i947i 

887i897i907i- 

91tt927i937i 

84o86o87o- 

91ol04ol05o 

53o54o60o- 

94a95o96o 

70o84o85o- 

97o98c99o 

4(90) 

887i897i907i- 

937i947i957i 

887i897i907i- 

91tt927i937i 

86o87o88o- 

91cl04cl05a 

55o57o60o- 

92o94o96o 

70o71o85o- 

97o98o99o 

5(86) 

847i857i867i- 

887i897i907i 

837i857i867i- 

877i887i897i 

82o83o84o- 

87o99ol00o 

50o52o54o- 

91o92o93o 

57o66o81o- 

93o94o95o 

6(86) 

837i857i867i- 

887i897i907i 

847i857i867i- 

877i887i897i 

80o82o83o- 

87o99ol00o 

50o52o54o- 

91o92o93o 

66o67o81o- 

93o94o95o 


NAAS Rating: 3.00 - Articles can be sent to editor @impactjournaIs.us 



Quantum Chemical Study of the Relationships Between Electronic Structure and Antiviral Activities Against 57 

Influenza a H1N1, Enterovirus 71 and Coxsackie B3 Viruses of Some Pyrazine-1,3-Thiazine Hybrid Analogues 


7(86) 


8 ( 86 ) 


9(86) 


10(86) 

11(93) 


12(93) 


13(93) 


14(90) 


15(90) 


16(90) 


17(86) 


18(86) 


19(86) 


837i847i857i- 

897t907t9l7t 

83jt857i867i- 

887t897t907t 

76a857i867i- 

887t897t907t 
847r857i867r- 
897t907i9l7t 
877i887i897i- 
947t977t987t 
887i897i937t- 
947t987t997t 
887i897i907i- 
947t957t967t 
877t897i907i- 
9l7i927t937t 
887i897i907i- 
9l7t927t937t 
877i887i897i- 
937t947t97a 
837i857i867i- 
887t897t907t 
847i857i867i- 
887t897t907t 
837i847i857t- 
897i:907r9 Itt 


847i857i867i- 

877r887r897i 

827i857i867i- 

877t887t897t 

847i857i867i- 

877t887t897t 

847i857i867i- 

877t887t897t 

9l7i927i937i- 

947t977t987t 

9l7i927i937i- 

947t957t967t 

907i9l7i937i- 

947t957t967t 

887i897i907i- 

9l7t927t937t 

877i897i907i- 

9l7t927i937t 

877i897i907i- 

9l7t927t937t 

847i857i867i- 

877t887t897t 

847i857i867i- 

877t887t897t 

837i857i867i- 

877i887i897i 


Table 10: Contd 

82g83g84g- 

87g99g101g 

80g82g84g- 

87g100g101o 

82o83g84g- 

87g100g101g 

80g83g84g- 

87o99g100g 

90g91g92g- 

95g107g108g 

90g91o92g- 

95g107g108g 

90g91g92g- 

95g107g108g 

84g87g88g- 

91g103g104g 

87o88g89g- 

91g103g104g 

84g87g88g- 

91g104g105g 

80o83g84g- 

87g99g100o 

82c83g84g- 

87g99g100g 

80o83g84g- 

87g99o100g 


52g55g56g- 

88g90g92g 

51g53o54g- 

91g92g93g 

52o55g57g- 

90g91g92g 

52g53g57g- 

88g90g92g 

57g59o60o- 


66g67g81g- 

93g94g95g 

61o67o81g- 

93g94g95g 

67o80g81g- 

93g94g95g 

67g69g81g- 

93o94o95g 

70c72g87o- 


96g99g100o 101g102g103g 

57c61g63o- 65o70c72c- 


98o99g100o 101g102g103g 

59o61g63o- 70o72c87c- 


95o99g100o 

56o58c61o- 

95g96g98g 

55g56g58g- 

95g96g98o 

58o60g70g- 

92g95g96g 

53g56g58g- 

91g92g93g 

52g53g57g- 

90g91g92g 

52g53g57g- 

90g91g92g 


101g102g103g 

70o84o85o- 

98o99g100g 

62o70g84g- 

98o99g100o 

66o70g86g- 

98o99g100o 

67o80g81o- 

93g94g95o 

67o80g81g- 

93o94g95g 

67g69g81g- 

93g94g95o 


DISCUSSIONS 

A Equations 1-3 show that there is a linear relationship between the variation of the biological property 
(log (IC 5 o)) and the variation of the numerical values of a definite set of local atomic reactivity indices. In the following we 
shall discuss case by case. 

Discussion of the results for the inhibition of the virus-induced cytopathic effect in Madin-Darby canine 
kidney cells infected with the H1N1 virus. 

Table 2 shows that the importance of variables in Eq. 1 is Sn N (LUMO+2)*» Fi 6 (LUMO+2)*> F 3 
(HOMO-2)*> Fi 8 (HOMO-2)*> Fi 0 (LUMO+2)*. Eq. 1 shows that a high inhibitory activity is associated with large 
numerical values for F 3 (HOMO-2)* and F 18 (HOMO-2)*, and with small numerical values for F 16 (LUMO+2)* and Fi 0 
(LUMO+2)*. As S 1 i N (LUMO+2)* is negative, then a high inhibitory activity is associated with large numerical values for 
this index. Atom 11 is a carbon in ring A (Figure 2). (LUMO) n*, (LUMO+l) u * and (LUMO+2) n * are 71 MOs (Table 9). 
Large negative values for Sn N (LUMO+2)* are obtained by shifting upwards the (LUMO+2) n * eigenvalue. Therefore, a 
good inhibitory activity is associated with empty MOs with low reactivity. This suggests that atom 1 1 is interacting as an 
electron donor. Atom 16 is a carbon in ring C (Figure 2). (LUMO+2) 16 *, (LUMO+l) 16 * and (LUMO)i 6 * have a 7t nature 
(Table 9). From Table 9 we can observe that (LUMO)ie* does not coincide with the molecule’s LUMO. As a low value for 
Fi 6 (LUMO+2)* is associated with a high inhibitory activity, a molecule with its five or six lowest empty MOs are not 
localized on atom 16 will be a good candidate. This suggests that atom 16 is interacting as an electron donor. Atom 3 is a 
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carbon in ring A (Figure 2). Large numerical values of F 3 (HOMO-2)* are associated with high inhibitory activity 
(let us remember that the values of the Fukui indices are in the [0,2] interval). Table 8 shows that (HOMO-2) 3 *, 
(HOMO-l) 3 * and (HOMO) 3 * have ji or o natures. This is a direct effect of the isopropyl substituent. The only way to 
integrate these facts is by suggesting that atom 3 could interact with two sites, one having it MOs and the other o MOs. 
Atom 18 is a carbon in ring C (Fig. 2). Large numerical values for Fi 8 (HOMO-2)* are associated with a high inhibitory 
activity. (HOMO-2 )i 8 *, (HOMO-1)i 8 * and (HOMO)is* have a jt nature (Table 10). This suggests that atom 18 is acting as 
an electron donor. Atom 10 is nitrogen in ring B (Figure 2). A high inhibitory activity is associated with small numerical 
values for F 10 (LUMO+2)*. Almost all (LUMO) 10 * have a n nature, all (LUMO+l)i 0 * have a ji nature and almost all 
(LUMO+l)io* have a 7t nature (Table 8). In quantum-chemical calculations of complex molecules like these ones, it is 
unavoidable that one or more empty MOs be localized on this atom. In this case we suggest the ideal situation is when the 
lowest empty local MOs do not correspond to the lowest empty MOs of the molecule. Regarding this point, it is suggested 
that quantum-chemical calculations be made with different substituents at position 12 (Figure 2) to examine their effects on 
the localization of the empty MOs on atom 10. Within the framework, it seems that atom 10 is acting as an electron donor. 
All the suggestions are displayed in the partial 2D pharmacophore of Figure 6. 



Figure 6: Partial 2D Pharmacophore for the Inhibition of the Virus-Induced Cytopathic Effect in 
Madin-Darby Canine Kidney Cells Infected with the H1N1 Virus 

Discussion of the Results for the Inhibition of Enterovirus 71 (EV71) 

Figure 4 shows that several points are outside the 95% confidence interval. A possible explanation is that an 
unknown factor has not been considered. Therefore, the following analysis must be taken with caution. Table 4 shows that 
the importance of variables in Eq. 2 is S 2 4 N (LUMO+2)*> Si 3 N (LUMO+2)*» Q l9 A high inhibitory activity is associated 
with a low negative net charge on atom 19. If S 2 4 N (LUMO+2)* is positive, a high inhibitory activity is associated with 
small numerical values for this index. If S 33 N (LUMO+2)* is positive, a high activity is associated with large numerical 
values. Atom 19 is the nitrogen of the NH 2 group attached to ring B. The requirement of a low negative net charge can be 
understood by suggesting that this atom is participating in a hydrogen bond of the form N-H...X, where X is an atom less 
electronegative that N or it is participating in an H-bond acceptor through its lone pairs. Atom 24 is one of the hydrogen 
atoms bonded to the abovementioned NH 2 group. Table 10 shows, as expected, that all occupied and empty local MOs are 
very far from the molecule’s frontier MOs. All they have a o nature (Table 10). Small positive numerical values for 
S 2 4 N (LUMO+2)* are associated with a high inhibitory activity. Low values are obtained by moving upwards the energy of 
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the associated eigenvalue, making this MO less reactive. The plot of S 24 N (LUMO)* vs. log (IC 50 ) (not shown) shows the 
same trend. The only explanation we have for these facts is that this H atom is facing a region with empty cr MOs and it is 
engaged in a repulsive interaction with them (zero electrons). Atom 13 is the carbon in ring C (Figure 2). A high inhibitory 
activity is associated with large positive values of S 13 N (LUMO+2)*. Table 9 shows that the three lowest empty local MOs 
have a n nature. Large positive values are obtained by shifting downwards the energy of the corresponding eigenvalue, 
making this MO more reactive. Therefore, we suggest that atom 13 seems to act as an electron acceptor. All the 
suggestions are displayed in the partial 2D pharmacophore of Figure 7. 



Figure 7: Partial 2D Pharmacophore for the Inhibition of Enterovirus 71 

We have stated that “it is important to stress that our hypothesis covers multi-step (for example, in the n-th step 
molecules must cross a pore) and multimechanistic (for example, to cross the pore molecules must interact consecutively 
with j unknown sites) processes. Therefore it seems logical to state that a necessary condition to obtain good 
structure-activity relationships is that all the steps and all the mechanisms inside each step must be the same for all the 
group of molecules under study ”. It seems that for this case the results suggest that all the molecules have not the same 
action mechanism. 

Discussion of the Results for the Inhibition of Coxsackievirus B3 (CVB3) 

Table 6 shows that the importance of variables in Eq. 3 is F 17 (LUMO+2)*> Si 8 E (HOMO-2)*> S 2 4 N 
(LUMO+l)*> Q 2 > Qi6> Q14. The analysis of Eq. 3 shows that a high inhibitory activity is associated with large values for 
F 17 (LUMO+2)*, large negative values for Si 8 E (HOMO-2)*, large positive values for S 24 N (LUMO+1)* 
(if the numerical values of this index are positive), negative values for Q 2 and Q 14 , and with positive values for Q 16 . Atom 
17 is a carbon in ring C (Figure 2). Table 10 shows that the three lowest empty local MOs have a ?t nature in all molecules. 
As a high inhibitory capacity is associated with large values for F [7 (LUMO+2)*, it is suggested that this atom is acting as 
an electron acceptor. Atom 18 is a carbon in ring C (Figure 2). Table 10 shows that the three highest occupied local MOs 
have a n nature. Large negative values for Si 8 E (HOMO-2)* are obtained by shifting upwards the energy of the 
corresponding eigenvalue, making this MO more reactive. This suggests that atom 18 is acting as an electron donor. 
Atom 24 is hydrogen (Figure 2). Table 10 shows that all MOs have a o nature. For positive values of S 24 N (LUMO+1)* a 
high inhibitory activity is associated with large values for this index. Large values are obtained by shifting downwards the 
energy of the associated eigenvalue, making this MO more reactive. This means that this hydrogen atom is participating in 
a H-bond of the N...H...X kind. Atom 2 is a carbon in ring A (Figure 2). A high inhibitory activity is associated with a 
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negative net charge. The data employed in LMRA shows that atom 2 has a net charge of about -0.4. Therefore it is 
suggested that atom 4 is engaged in an electrostatic interaction with a counterpart having a positive net charge. Atom 16 is 
a carbon in ring C (Figure 2). A high inhibitory capacity is associated with positive values for its net charge. The analysis 
of the numerical values in the data matrix shows that the net charge of this atom is close to zero (-O.Ox) in all molecules. 
Therefore a substitution producing positive net charges should help to enhance activity. This suggests that atom 16 is 
facing a negatively charged group. Atom 14 is a carbon in ring C (Figure 2). A high inhibitory activity is associated with 
negative values for its net charge. Almost net charges are negative in the data matrix employed for LRMA 
(one net charge is positive). Therefore we suggest that atom 14 is facing a positively charged group. All the suggestions are 
displayed in the partial 2D pharmacophore of Figure 8. 



Figure 8: Partial 2D Pharmacophore for the Inhibition Coxsackievirus B3 


In quantum chemical calculations many times some eigenvalues of empty MOs are negative. This is the reason 


why in our analysis we always stated that “/or a positive value of Si (MO) a high activity is associated with”, etc. 


The reader will easily verify that for negative values the analysis is exactly the same. 


CONCLUSIONS 


In conclusion, we have obtained at least two statistically significant equations relating the variation of the 
inhibitory activity against coxsackievirus B3 and H1N1 virus with the variation of the values of definite local atomic 
indices of a common skeleton. The results for Enterovirus 71 suggest that all the molecules studied here do not seem to 
have the same mechanism of action against this pathogen. 
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